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ABSTRACT 

We report observations of HD 80606 using the 10.4-m Gran Telescopio Canarias (GTC) 
and the OSIRIS tunable filter imager. We acquired very-high-precision, narrow-band 
photometry in four bandpasses around the K I absorption feature during the January 
2010 transit of HD 80606b and during out-of-transit observations conducted in January 
and April of 2010. We obtained differential photometric precisions of ~ 2.08 x 10 -4 
for the in-transit flux ratio measured at 769.91-nm, which probes the K I line core. 
We find no significant difference in the in-transit flux ratio between observations at 
768.76 and 769.91 nm. Yet, we find a difference of ~ 8.09 ±2.88 x 10~ 4 between these 
observations and observations at a longer wavelength that probes the K I wing (777.36 
nm). While the presence of red noise in the transit data has a non- negligible effect on 
the uncertainties in the flux ratio, the 777.36-769.91 nm colour during transit shows 
no effects from red noise and also indicates a significant colour change, with a mean 
value of ~ 8.99 ±0.62 x 10 -4 . This large change in the colour is equivalent to a ~ 4.2% 
change in the apparent planetary radius with wavelength, which is much larger than 
the atmospheric scale height. This implies the observations probed the atmosphere 
at very low pressures as well as a dramatic change in the pressure at which the slant 
optical depth reaches unity between ~770 and 777 nm. We hypothesize that the excess 
absorption may be due to K I in a high-speed wind being driven from the exoplanet's 
exosphere. We discuss the viability of this and alternative interpretations, including 
stellar limb darkening, starspots, and effects from Earth's atmosphere. We strongly 
encourage follow-up observations of HD 80606b to confirm the signal measured here. 
Finally, we discuss the future prospects for exoplanet characterization using tunable 
filter spectrophotometry. 

Key words: planetary systems - stars: individual (HD 80606) - techniques: photo- 
metric 



1 INTRODUCTION 

Discoveries of extrasolar planets which transit their host 
star provide valuable opportunities to measure the physical 
properties of exoplanetary atmospheres. The physical char- 
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acteristics of an exoplanetary atmosphere can be probed by 
transmission s pectroscopy observed aga i nst the spectru m of 
the host star. ISeager fc Sasselovl (^POOp . iBrownl i|200ll ) and 
iHubbard et al.l (|200ll ) developed models that predicted such 
absorption, particularly from Na I, K I, and other alkali 
metals. Subsequent refinements of such models have con- 
firmed that in the optical wavelength regime, the strongest 
lines are expected from the Na I resonance lines (AA 589.6, 
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589.0 nm) and the K I resonance lines (AA 769.9, 766.5 nm) 
(e.g., iBarmanl [20071 ; iFortnev et all l201Cf) Fl In the optical, 
the cores of the atomic features of Na I and K I are rela- 
tively narrow. For this reason, medium to high resolution 
spectrographs can be used to compare the in-transit stellar 
spectrum to the out-of-transit stellar spectrum. The absorp- 
tion of stellar photons in the exoplanetary atmosphere leads 
to excess absorption in the in-transit stellar spectrum when 
compared to the out-of-transit spectrum. In photometric ob- 
servations, this leads then to deeper transits and a larger 
apparent size of the planet at the absorbing wavelengths 
llBrownl | 200 If), with v ariations of order the atmospheric scale 
height l|Fortnevll2005l ). Such measurements in strong optical 
transitions can also constrain the atmospheric metallicity, 
rainout of condensates, distribution of absorbed stellar flux, 
and photoionization of atmospheric constituents. 

The first detection of absorption due to an exo- 
planetary atmosphere came from Na I observations of 
HD 209458b using the Space Telescope Imaging Spectro- 
graph (STIS) on board th e Hubble Space Telescope (HST) 
i|Charbonneau et all [2002] ) . Unfortunately, the subsequent 
failure of the STIS instrument prevented similar observa- 
tions for more than 5 years. Thus, attention was directed 
toward making such observations from the ground (e.g . , 
iMoutou et aLI l200ll ; IWinn et all 12004 iNarita et all 120051) . 
The second detection of absorption due to an exoplanetary 
atmosphere, this time from the ground, was also made of 
Na 1 in observations of H D 189733b using the 9.2-m Hobby- 
Eberly Telescope (HET) jRedneld et alj|200Sl) . Further de- 



tections of Na 1 in the atmosphere of HD 209458b were 



((Snellen et al.l 


20081). from HST bv 


Sing et all l(2008a!) and 


from Keck by 


Langland-Shula et al. 


(2009). The recent re- 



pair of STIS and installation of the Cosmic Origins Spec- 
trograph (COS) on board HST has enabled new optical 
and ultraviolet (UV) transmission spectrum observations 
of exoplanetary atm ospheres, extended exospheres, and au- 
roral emission (e.g ., ILinskv et aill2010l ; iFossati et aill2010l ; 
iFrance et aHl2010h . 

Compari ng the surprisingly weak Na I absorption in 
HD 209458b i|Charbonneau et al.l2002j ; lKnutson et al.ll2007l ) 
to the 3 times stro nger Na I absorption of HD 189733b 
l|Redfield et al.|[2008l ) suggests that the two planets have dif- 
ferent atmospheric structures. Theorists have suggested nu- 
merous mechanisms such as adjustments to the metallicity, 
rainout of condensates, distrib ution of absorbed stellar flux, 
or ph otoionization of sodium (lBarman| 2007; Fo rtnev et al.l 
|2003| ). In particular, iBarman et alj ( 20021 ) suggested that 
non-LTE Na level populations were the cause of the weak 
Na feature observ e d in HD 2094 58b, and a reana lysis of 
the lKnutson et al.l (|2007l ) data bv lSing et ail l|2008al fbT) sug- 
gested that Na condensation or Na photoionization in HD 
209458b's atmosphere was the best explanation for match- 
ing the data, given the Na line shapes they derived. It is 



1 We caution that these lines are most prominent for hot-Jupiter- 
like planets with a certain range of atmospheric temperatures. 
Atmosphere models generated for H D 80606b at the time of tran- 
sit [based on IFortnev et"al] <20ld) l do not predict a significant 
K I absorption feature, due to the low equilibrium temperature 
of 500 K. We refer the reader to i|4,4l for further discussion. 



clear that comparisons of the atmospheric properties of dif- 
ferent transiting planets will be critical to understanding the 
atmospheric properties of exoplanets as a whole. Although 
still small, the list of detected atoms and molecules is grow- 
ing. In addition to Na I, several molecules have been de- 
tected, primarily in the infrared, with both space-based and 
ground-based platf o rms, i ncluding CO, CO2, H 2O, and CH4 
i|Swain et all I2008L 120091 ; ISnellen et alll2010j ). Other HST 
observations using the Advanced C amera for Surveys (ACS) 
did not detect K I in HD 189733b l|Pont et al.ll200Sl ). If de- 
tections of constituents in the extended exosphere are in- 
cluded, then H 1, C 11, O 1, Mg 11, and other metals have also 
been detected (IVidal-Madiar et alj|2003 . |2004 ILinskv et ail 
|2010| ; IFossati et alj|2010h . 

Each new detection provides not only compositional in- 
formation, but another window into the physical properties 
of the exoplanetary atmosphere (e.g., condensation, wind 
speed, photoionization). Even though atmosphere models do 
not predict a significant K I feature in HD 80606b, it remains 
of great interest to observationally determine the level of K I 
absorption in its atmosphere, since K I is generally predicted 
to be the second strongest transmission spectrum signature 
in the optical wavelength range. Further, Na I and K I probe 
different layers of the atmosphere. Measurements of K I can 
test the hypothesis that the low abundance of Na I on HD 
209458b may be due to a high altitude layer of clouds or 
haze. Finding low abundance for both Na I and K I would 
be consistent with either the cloud hypothesis or with the 
photoionization hypothesis, as both are very easy to ionize. 
Finding that only Na I is significantly depleted would point 
to alternative models with complex atmospheric chemistry 
(e.g., incorporation into grains, odd temperature structure, 
unexpected mixing patterns) . Finally, in principle, future ob- 
servations could probe temporal variability of Na I and K I 
due to high-speed, high-altitu de winds and/or diff erences in 
the leading and trailing limb IFortnev et al.ll2010l ). 

All of the above atmospheric studies were based on 
observations using high-resolution spectrographs. Here, we 
describe a new technique that utilizes fast, narrow-band, 
spectrophotometry with the Optical System for Imaging 
and low Resolution Integrated Spectroscopy (OSIRIS) in- 
stalled on the 10.4-m Gran Telescopio Canarias (GTC) to 
probe the composition and other properties of the atmo- 
spheres of exoplanets that transit bright stars (see Ej2j) . Fast 
line spectrophotometry can be much more efficient (e.g., 
~ 34% with GTC/OSIRIS) than typical high resolution 
spectrographs (~l-2%) thanks to the use of a tunable fil- 
ter (TF) rather than diffraction gratings. Further, this tech- 
nique has the potential to be less sensitive to several sys- 
tematic noise sources, such as seeing variations that cause 
line- variations in wide spectrograph slits (specifically in non- 
fiber fed spectrographs), atmospheric variations (since ref- 
erence stars will be observed simultaneously), and/or fiat- 
fielding errors (since on- and off-line data are obtained at 
the same detector location). Thus, spectrophotometry with 
a TF technique is particularly well-suited for observing a 
narrow spectral range of atomic absorption features, with- 
out suffering from the inefficiencies or potential systematic 
uncertainties of high-resolution spectrographs. 

Here we present results of such observations of the Jan- 
uary 2010 transit of HD 80606b using the GTC and the 
OSIRIS TF imager. HD 80606b was originally discovered 
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by radial velocity observations (|Naef et al 1l200lh and was re- 
markable due to its very high eccentricity (e = 0.93). Only 
several years later did Spitzer and ground-based observa- 
tions reveal that t he planet passes both behind and in front 



of its host star dLaughlin et alj 120091 ; Fossev et all 20091 ; 



iGarcia-Melendo fc McCulloughl 120091 ; iMoutou et al J 120091 ) 

Spectroscopic observations revealed that the angular mo- 
mentum axis of the stellar rotation and that of the or- 



bital planet are misalig ned (|Moutou et al. I l2009l ; iPont et al.l 
2009; IWinn et alj 12009 ). Given the infrequent transits and 
long transit duratio n (~ 12 hours), fol l ow-up observations 
are quite challenging. IWinn et al.l (|2009r i. lHidas et al.l (|2010l 'l 
and IShporer et all (|20ld ) were able to characterize transits 
of HD 80606b with longitudinally distributed networ ks of 
ground-based observatories, and Iflebrard et al.l (|2010f ) ob- 
served the January 2010 transit using the Spitzer spacecraft. 

The Spitzer observations constr ain the thermal prop- 
erties of the plane t's atmosphere l|Laughlin et al.l 2009; 
iHebrard et afll201Ch . To the best of our knowledge, the ob- 
servations presented here are the first to attempt to detect 
atmospheric absorption by HD 80606b. While existing atmo- 
sphere models predict that HD 80606b would not have any 
significant K I feature due to its high surface gravity and 
cold atmosphere at the time of transit (e.g., see H4.4|) . our 
observations test this prediction. Even though models do not 
predict a K I feature, exoplanet observations have a track 
record of unexpected discoveries. Furthermore, in principle, 
depending on the atoms/molecules found in the atmosphere, 
these observations could yield information about how the 
planet cools, independent of any observations of the ther- 
mal phase curve of this system. In principle, transmission 
spectroscopy also provides a way to characterize transiting 
planets in eccentric orbits, which either do not pass behind 
their host star or which are too cool to detect via occultation 
when they do pass behind the star. 

Finally, we note that HD 80606 is one of the best sys- 
tems for making very precise spectrophotometric measure- 
ments. HD 80606 is the brightest of the transiting planet 
host stars which have a comparably bright reference star 
very nearby (~20 arcsec). Also, the long duration between 
the 2nd and 3rd points of contact (~6 hours) of HD 80606b 
provides time to collect a large amount of in-transit data in 
a single transit. Thus, we expect that all else (e.g., observing 
conditions) being equal, HD 80606b permits the most pre- 
cise spectrophotometric measurements of any known system 
(at least with observations of a single transit). 

This paper presents extremely precise measurements of 
the variation in HD 80606b's apparent radius with wave- 
length near the K I feature, which in turn can help us test 
the predictions of atmosphere models. Sj2] describes our ob- 
servations and data analysis procedures. We describe the 
results of our observations in fJ3] in 3Uwe interpret the re- 
sults, and we summarize our conclusions and discuss the 
future prospects for the method in fj5]and Sj6] 



2 OBSERVATIONS 

HD 80606 and its nearby companion (HD 80607) are both 
bright G5 dwarves of a similar magnitude (V~9) and colour. 
On three nights, we measured the flux of both HD 80606 
(target) and HD 80607 (reference) simultaneously. We cy- 



cled through a set of four wavelengths throughout the ob- 
servations. On the night of 2010 January 13-14, the planet 
was in transit for the duration of our observations, and we 
measure an "in-transit" flux ratio of HD 80606 to HD 80607 
for each wavelength. We repeated the observations on 2010 
January 15 and 2010 April 4, when the planet was not tran- 
siting HD 80606, allowing us to measure the "out-of-transit" 
(OOT) flux ratio of HD 80606 to HD 80607 for each wave- 
length. Our results (fJ3| are based on the ratio of in-transit 
flux ratio (target over reference) to out-of-transit flux ratio 
(target over reference). Any changes in the Earth's atmo- 
sphere from one night to the next should affect both the 
target and reference star similarly. By making differential 
measurements of the colour during the same transit and at 
similar atmospheric conditions, this method allows for ex- 
tremely precise measurements of the transit depth at dif- 
ferent wavelengths. While night-to-night variability in the 
atmospheric conditions or either of the stars could cause a 
systematic scaling of the transit depth measurements, the 
relative wavelength dependence of the apparent planet ra- 
dius is largely insensitive to either of these potential sys- 
tematics. We refer the reader to i]4.9.1l and £|4.9,3l for further 
discussion. 



2.1 In- Transit and Out-of- Transit Observations 

We observed a partial transit of HD 80606b on 2010 Jan- 
uary 13-14 and acquired baseline data on 2010 January 15 
and 2010 April 4 to establish the OOT flux ratios. For our 
observations, we used the TF imaging mode of the OSIRIS 
instrument installed on the 10.4-m GTC, which is located at 
the Observat orio del Rogue de los M uchachos on the island 
of La Palma (|Cepa et alj|2000l . \200j ). In the TF mode, the 
user can specify custom bandpasses with a central wave- 
length of 651-934.5 nm and a FWHM of 1.2-2.0 nm. The 
effective wavelength decreases radially outward from the op- 
tical centre; because of this effect, we positioned the target 
and its reference star at the same distance from the opti- 
cal centre and on the same CCD chip. The observed wave- 
lengths described below refer to the location of the target 
(and reference) on the CCD chip. 

During the transit observations and baseline observa- 
tions on 2010 January 15, exposures of the target and its 
reference star cycled through four different wavelengths (all 
with a FWHM of 1.2-nm): one on the predicted core of the 
K I line (769.75-nm), one to the blue side (768.60-nm) and 
two redward of the K I feature (773.50 and 777.20 nm). As 
the tunings for the TF are set by the order sorter (OS) fil- 
ter used, our bluest wavelength is the bluest wavelength we 
could observe at in the wing of the K I line and still observe 
within the same OS filter as the "on-line" wavelength (i.e., 
at the location of the core of the K I line). We then chose 
two wavelengths redward of the K I line in order to sam- 
ple more of the structure/ wings around the K I line. The 
reddest bandpass was chosen since we expect to see (for a 
typical hot- Jupiter) a maximum difference between the flux 
ratio in the on-line bandpass and around that reddest band- 
pass. In order to maximize the signal-to-noise ratios in the 
on-line wavelength and in the reddest off-line wavelength, in 
each sequence we observed on-line three times, at the red- 
dest off-line wavelength two times and at the other off-line 
wavelengths one time each. During the transit, the observ- 
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ing sequence from the GTC was as follows: 769.75, 768.60, 

769.75, 773.50, 769.75, 777.20, 777.20 nm (repeat). 

We emphasize that these wavelengths were chosen to be 
around the location of the K I feature in HD 80606b's atmo- 
sphere. In order to observe on the K I feature (which has a 
rest wavelength of ~769.9 nm) in the frame of the planet, we 
accounted for the Doppler shifts due to the Earth's motion 
around the Sun, the system's radial velocity and the planet's 
non-zero radial ve locity during trans it [—59.6 km s _1 based 
on velocities from IWinn et al](|2009h ]. After accounting for 
these effects, the observed wavelengths in the frame of the 
planet are redshifted by 0.16 nm to 769.91 nm (on-line) and 

768.76, 773.66 and 777.36 nm (off-line). The observed wave- 
lengths in the frame of the star are essentially the same as 
observed on Earth due to the small systemic velocity of the 
HD 80606 planetary system and the Earth's small barycen- 
tric velocity on the night of the transit. For the remainder 
of the paper, we report the wavelengths as observed in the 
frame of the planet when discussing results from the transit 
observations. 

A similar sequence as described above was used for the 
baseline observations taken on 2010 April 4, but the ob- 
served wavelengths were corrected for the Doppler shift due 
to the planet's orbital velocity on that specific date (~23.9 
km s _1 ) in order to match the wavelengths observed during 
the transit. Thus, the wavelengths observed on 2010 April 4 
(from the GTC) are 770.00 nm (on-line) and 768.86, 773.76 
and 777.45 nm (off-line). 

Transit observations of HD 80606b began at 22:28 UT 
on 2010 January 13 (during ingress) and ended at 7:15 UT 
on 2010 January 14 (around the beginning of egress and 
including astronomical twilight), during which the airmass 
ranged from ~1. 08-1. 72. The observing conditions were pho- 
tometric, with a clear sky and a dark moon. No data was 
taken between 5:20 and 5:50 UT on 2010 January 14 due to 
recalibration of the TF during that time. The actual seeing 
varied between 0.7 and 0.9 arcsec during the transit obser- 
vations, but we used a slight defocus to increase efficiency 
and reduce the impact of pixel-to-pixel sensitivity variations. 
Therefore, the defocussed FWHM of the target varied from 
~0.9-2.3 arcsec (7-18 pixels) during the transit. For the por- 
tion of the light curve used in our analysis (see H2.2I) . the 
FWHM was much more stable than is indicated by the range 
given above, with a typical value between 10 and 14 pixels 
and a mean value of 12 pixels. Even with an autoguiding 
system, the target's centroid coordinates shifted by ~9-10 
pixels over the course of the night. We used lxl binning 
and a fast readout mode (500 kHz) to readout a single win- 
dow of 300 x 600 pixels (located on one CCD chip) in order 
to reduce the dead time between exposures. This window is 
equivalent to a field-of-view of ~ 38x76 arcsec, so the only 
stars in our field were HD 80606 and a single reference star, 
HD 80607. Each individual observation was followed by an 
average dead time of less than 4-s for readout and to switch 
between TF tunings. We used 10-s exposures, resulting in 
an overall cadence of about 14-s for each observation. Due 
to the short exposure time used, the sky background level 
was low enough that we did not need to discard any images 
taken during astronomical twilight. 

Baseline observations were taken from 5:50 to 7:10 UT 
(i.e., also through the beginning of astronomical twilight) on 
2010 January 15, but the data was highly scattered, so we do 



not include it in our primary analysis^ Additional baseline 
observations took place on 2010 April 4 from 21:30 (includ- 
ing the end of astronomical twilight) to 0:00 UT. The ob- 
serving conditions were photometric and taken during grey 
time, using the same set up as the in-transit observations de- 
scribed above. During the observations, the airmass ranged 
from ~1. 08-1. 20, and the actual seeing varied between 1.4 
and 1.6 arcsec (11-12.5 pixels), so the telescope was not 
intentionally defocussed. The target's centroid coordinates 
shifted by ~5-8 pixels during the observations. The expo- 
sure time was changed from the initial exposure time of 10-s 
to 8-s and then again to 11-s to counteract variations in the 
seeing as well as increasing airmass while avoiding satura- 
tion and maintaing a high number of counts. In our analysis, 
we discard the 10-s data because a majority of the images 
were saturated. We tested using the OOT flux ratios from 
the 8 and 11 s data individually in our analysis and found 
that they produced very similar results. Thus, we combine 
the 8 and 11 s data to establish the final OOT flux ratios 
(see i]2.2[) and to achieve the longest usable baseline possible. 

2.2 Data Reduction and Analysis 

Observations taken with OSIRIS prior to mid-March 2010 
suffered from a higher than expected level of dark current 
despite the short exposure times used. Therefore, we used 
standard IRAF procedures for bias and dark subtraction 
and flat field correction for the January 2010 transit obser- 
vations of HD 80606. We note that the flat fields for these 
observations did not produce the pattern of having the to- 
tal number of co unts in the d o me fla t fields decreasing with 
time as seen by I Colon et al.l l|2010h . so we use almost all 
(65 out of 75) dome flats for each observed wavelength in 
our analysis (the 10 dome flats not included in the analysis 
were overexposed). A new dewar fixed the problems with the 
dark current before the April 2010 observations took place, 
so for the baseline data we performed standard bias subtrac- 
tion and flat field correction (combining all 133 flats taken 
for each observed wavelength) and did not need to subtract 
dark frames. 

Because of the very small readout window used for 
our observations, our images do not contain the sky (OH) 
emission rings that occur due to the TF's small band- 
pass and position-dependent wavelength. Therefore, we per- 
formed simple aperture photometry on the target and refer- 
ence star using the standard IDL routine APER[j for a range 
of aperture radii. We measured the root mean square (RMS) 
scatter of the flux ratio (equal to the target flux divided 
by the reference flux) at the bottom of the transit (for the 
January 2010 data) and for the individual 8 and 11 s data 
taken OOT (in April 2010) in each bandpass. We consid- 
ered the results for each bandpass and adopted an aperture 
radius of 28 pixels (3.6 arcsec) for the in-transit data and 
32 pixels (4.1 arcsec) for the OOT data, as these were the 
aperture radii that typically yielded the lowest RMS scatter. 
The radii of the sky annulus used for the reduction of both 
data sets was 68 to 74 pixels in order to completely avoid 
any flux from the target or reference star. We have included 

2 See i|4.2l for further discussion of this data set. 

3 http://idlastro.gsfc.nasa.gov/ 



Probing the atmosphere of HD 80606b 5 



Table 1. Absolute Transit Photometry from 2010 January 13 



A (nm) 



HJD 



^target ^*re/ 



768.76 2455210.4428 2789803 2501667 



769.91 2455210.4414 1207952 1081769 



773.66 2455210.4419 1789988 1603033 



777.36 2455210.4423 2441876 2188980 



The wavelengths included in the table are the observed wave- 
lengths in the frame of the planet (see text for additional details). 
The time stamps included here are for the times at mid-exposure. 
Ftarget and F re f are the absolute flux measurements of HD 80606 
and HD 80607. The full table is included online, while a portion 
is shown here so the reader can see the formatting of the tabic. 



Table 2. Absolute OOT Photometry from 2010 April 4 



A (nm) t exp (s) 



HJD 



Ftar 



768.1 



2455291.4245 5550569 4951634 



770.00 8 2455291.4243 5665554 5060730 




Time from Mid-Transit (Hours) 



Time from Mid— Transit (Hours) 



Figure 1. Absolute fluxes of HD 80606 (a) and HD 80607 (b) as 
measured on 2010 January 13-14. The different light curves repre- 
sent the fluxes as measured nearly simultaneously in the different 
bandpasses, with the black, blue, brown and red light curves rep- 
resenting the 769.91, 768.76, 773.66 and 777.36 nm data. These 
data have not been corrected for airmass or decorrelated in any 
way. Note the break in the data around two hours after mid- 
transit due to recalibration of the tunable filter. The vertical solid 
lines indicate the expected beginning and end of the transit, and 
the vertical dotted lines indicate the end of ingre ss and the be- 
ginnin g of egress [based on durations estimated bv | Hebrard et al.l 
]|201Ch and the transit ephemeris from lShporer et alj l|201ol 'l]. The 
vertical dashed lines indicate the ~4 hour interval around mid- 
transit that our analysis focused on (see text for further details). 



773.76 



2455291.4249 6125160 5466145 



777.45 



2455291.4253 7006625 6252784 



Columns are similar to Table [JJ except the wavelengths included 
in the table are the wavelengths as observed from the GTC (see 
text for additional details). The second column contains the ex- 
posure time for the observations, as observations based on two 
different exposure times were included in our analysis. The full 
tabic is available online, and a portion is shown here so the reader 
can see the formatting of the table. 



the results of our aperture photometry in Tables [TJ and [2] 
and illustrate the results in Figures [TJ and [2] As illustrated, 
the flux in each bandpass displayed large variations during 
parts of the observations (particularly during parts of the 
transit), and we take this into consideration in our analysis 
(see £TTJ|. 

We present the raw in-transit light curves in Figure 
[3] which were computed by dividing the flux in the target 
aperture by the flux in the reference star aperture and then 
normalizing by the weighted mean OOT flux ratio for each 
bandpass (see Sj3]for details on the computation of the mean 
flux ratios). In an attempt to reduce systematic trends seen 
in our transit light curves, we applied the ext ernal param- 
eter decorr elation (EPD) technique (see, e.g., iBakos et alj 
120071 l20ld ) to each transit and baseline light curve. Note 
that for the transit light curve, we only applied EPD to the 
~4 hours centered around mid-transit, or 3:36 UT on 2010 
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Figure 2. Similar to Figure [TJ but for the out-of-transit data 
taken the night of 2010 April 4. Note that the discontinuity in 
the fluxes around 2455291.48 is due to a change in the exposure 
time (from 8-seconds to 11-seconds). 



January 14, as estimated bv lShporer et al. (2010)Q Specifi- 
cally, we decorrelated each individual light curve against the 
following parameters: the centroid coordinates of both the 



4 T his ephemeris i s in be tween that given bv lWinn et al. (2009) 
and lHebrard et al ] [|2010h . The choice of ephemeris used does not 
significantly affect our results. 
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Figure 3. Transit light curves as observed nearly simultaneously 
in different bandpasscs on 2010 January 13-14. The on-line light 
curve (769.91-nm) is shown in black, and the off-line light curves 
(768.76, 773.66 and 777.36 nm) are shown in blue, brown and 
red. The flux ratio for each bandpass has been normalized to the 
weighted mean OOT flux ratio estimated from the baseline data 
acquired in April 2010, but the data has not been corrected for 
airmass or dccorrelated in any way. The off-line light curves have 
been arbitrarily offset by 0.006, 0.012 and 0.018, and error bars 
are not shown for clarity. The vertical solid, dotted and dashed 
lines are the same as in Figure [l] 

target and reference, the sharpness of the target and ref- 
erence profiles [equivalent to (2.35/FWHM) 2 ] and the air- 
mass. As illustrated in Figure 21 EPD removed most of the 
correlations in the in-transit data. For reference, we show 
the correlations between the in-transit data and the target's 
FWHM and centroid coordinates both before and after EPD 
has been applied in Figure [5] For the baseline data, we per- 
formed EPD for the 8 and lis data series separately, but we 
then combined the two data sets to compute the weighted 
mean flux ratio and its uncertainty for each bandpass as de- 
scribed in fj3] The results of the decorrelation for the OOT 
data are illustrated in Figure[6] As a result of applying EPD, 
the RMS scatter in each of the bandpasses improved by as 
much as ~25%, but decorrelating the light curves against 
the above parameters did not completely remove the sys- 
tematics that are seen in our data. In a further attempt to 
remove systematics, we also tried a quadratic decorrelation 
against the sharpness of the target and reference profiles, as 
that was the only parameter that showed a possible resid- 
ual systematic pattern after EPD was applied. However, the 
quadratic decorrelation did not reduce systematics in our 
light curves any further. We discuss other potential sources 
of systematics in detail in &14.9I 

Because our goal is to compare the depths of the transit 
in each bandpass, the rest of our analysis focuses on the data 
from the bottom of the transit as presented in Figure [3] and 
highlighted in Figure [7] i.e., the ~4 hours centered around 
mid-transit. Note that the light curves shown in Figure [7] 
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Figure 4. Relative in-transit flux ratio normalized to the relative 
OOT flux ratio as measured on 2010 April 4. The relative flux 
before (a) and after (b) EPD was applied is shown. The different 
colors represent the flux ratios as measured in the different band- 
passes, with the colors the same as in Figure [3] Note that EPD 
was only applied to the ~4 hours centered around mid-transit 
(i.e. the bottom of the transit light curve). The data shown has 
not been binned, but the different light curves have been offset 
arbitrarily for clarity. 

have been corrected using EPD. We also discarded points 
that had a flux ratio greater than 3<r from the mean of the 
bottom of the transit light curve. This resulted in discarding 
four points from the reddest light curve (777.36-nm). We also 
discarded several exposures from each wavelength that were 
unusable due to saturation. The different panels in Figure 
[7] illustrate the deviation between the magnitude of the on- 
line flux ratios and each of the off-line flux ratios, which will 
be discussed in detail in Sj3]and 

We estimated the uncertainties in the flux ratios by 
computing the quadrature sum of the photon noise for HD 
80606 and HD 80607, the uncertainty in the sum of the sky 
background (and dark current, for the in-transit observa- 
tions) and the scintillation noise for the two stars. We as- 
sume Poisson statistics to compute the uncertainty in the 
sky background, and the noise due to scintillation was es - 
timated from the rela tion given by iDravins et all (|l998h . 
based on lYound (|l967l ). We caution that this empirical re- 
lation might overestimate scintillation for large telescopes 
located at excellent sites such as La Palma. Regardless, the 
relation demonstrates that scintillation is still a small con- 
tribution to the total error budget for these observations. 
The flat field noise is also negligible compared to the pho- 
ton noise, so we do not include it in our determination of 
th e measu r ement uncertainties. Based on the relation given 
by iHowelll (|2006t ). which computes the standard deviation 
of a single measurement in magnitudes and includes a cor- 
rection term between the error in flux units and the error 
in magnitudes, we find the median total uncertainties in 
the flux ratio for each exposure to be 0.538, 0.532, 0.514 
and 0.486 mmag at 768.76, 769.91, 773.66 and 777.36 nm 
(over the bottom of the transit). The RMS of the transit 
light curve is comparable, but slightly larger, with values of 
0.585, 0.667, 0.631 and 0.662 mmag for those wavelengths. 
The median total uncertainties for the OOT observations 
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Figure 5. Correlations between the normalized in-transit flux 
ratio and the target FWHM and x and y centroid coordinates, 
before (lefthand column) and after (righthand column) EPD has 
been applied. All four bandpasses are shown in each panel, with 
the colours the same as in Figure [3] Similar results were obtained 
when decorrelating the data against the reference parameters but 
are not shown here. 
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Figure 6. Relative OOT flux ratio as measured on 2010 April 
4. The relative flux before (a) and after (b) EPD was applied is 
shown. The different colors represent the flux ratios as measured 
in the different bandpasses, with the colors the same as in Figure 
\3\ Note the small break in the data around 2455291.48 where the 
exposure time was changed. The data has not been binned, but 
the different light curves have been offset arbitrarily for clarity. 



Figure 7. Corrected light curves for observations of the bot- 
tom of the transit as observed nearly simultaneously in differ- 
ent bandpasses on 2010 January 13-14. In each panel, the black 
points illustrate the measurements taken in the on-line (769.91- 
nm) bandpass. We also show measurements taken in each of the 
off-line bandpasses (768.76, 773.66, 777.36 nm) in each of the re- 
spective panels (a, b, c) for comparison to the on-line flux ratios. 
The data shown here has been decorrelated. The colours and nor- 
malizations are the same as in Figure [3] but no offsets have been 
applied. Here, we have binned the data and error bars simply for 
clarity. 



are calculated in a similar way, but the uncertainties for the 
8 and lis data sets were scaled by the flux ratios for each 
respective set in order to compute a weighted uncertainty. 
Thus, the median total (weighted) uncertainties in the flux 
ratio are 0.657, 0.650, 0.627 and 0.592 mmag while the esti- 
mated RMS is quite comparable, with values of 0.562, 0.605, 
0.554 and 0.586 mmag for 768.73, 769.87, 773.63 and 777.32 
nm. 

The complete photometric time series for each bandpass 
of the in-transit data (uncorrected and unnormalized) is re- 
ported in Table [3] while the photometric time series (both 
before and after EPD was applied) for the transit bottom 
and the April observations are reported in Tables [4] and [5] 
The weighted mean flux ratios for both the in-transit and 
OOT data (see ij3] for more details) are given in Table |6j 
along with their uncertainties. 
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Table 5. Relative OOT Photometry from 2010 April 4 



A (nm) t exp (s) HJD F raUo F ratio Uncertainty 

(raw) (corrected) 

768.86 8 2455291.4245 1.12096 1.12050 0.00058 



770.00 8 2455291.4243 1.11951 1.11888 0.00058 



773.76 8 2455291.4249 1.12056 1.11973 0.00056 



777.45 8 2455291.4253 1.12056 1.11983 0.00053 



The wavelengths included in the table are the wavelengths as observed from the GTC (see text for additional 
details). The time stamps included here are for the times at mid-exposure. The flux ratios are presented 
both before (raw) and after (corrected) EPD was applied. The full table is available online, and a portion is 
shown here so the reader can see the formatting of the table. 



Table 3. Relative Transit Photometry 



A (nm) HJD F ra tio Uncertainty 

768.76 2455210.4428 1.11518 0.00086 

769.91 2455210.4414 1.11665 0.00124 

773.66 2455210.4419 1.11663 0.00104 

777.36 2455210.4423 1.11553 0.00091 



The wavelengths included in the table are the observed wave- 
lengths in the frame of the planet (see text for additional details). 
The time stamps included here are for the times at mid-exposure. 
Fratio represents the relative flux ratio between the target and ref- 
erence star (i.e., Ft ar get / F re f). The full table is available online, 
and a portion is shown here so the reader can see the formatting 
of the table. 

3 RESULTS 

As illustrated in Figure [7] we can see by eye a hint of a 
deviation between the in-transit flux ratios observed at the 
on-line wavelength and the red off-line wavelengths, but no 
clear deviation is seen when compared to the bluest off-line 
wavelength. Despite evidence of time-correlated systemat- 
ics in our data, we emphasize that the error bars shown 
in Figure [7] are binned error bars, which illustrate that our 
measurement uncertainties are larger than any residual sys- 
tematics present in the light curves and that the deviations 
in the flux ratios between the different bandpasses are real. 
We refer the reader to our discussion of possible systematic 
sources in ^4.91 

In Figure [8] we plot histograms of the (unbinned) flux 
ratios at the bottom of each of the transit light curves, 
where the flux ratios have been normalized against the mean 
OOT flux ratio for each respective wavelength. These his- 
tograms further illustrate that the flux ratios for the on-line 



Table 4. Normalized Photometry from around Mid- Transit 



A (nm) HJD F ra tio F ra tio Uncertainty 

(raw) (corrected) 

768.76 2455210.5680 0.99033 0.98962 0.00054 

769.91 2455210.5670 0.99041 0.98912 0.00054 

773.66 2455210.5671 0.99036 0.98909 0.00052 

777.36 2455210.5675 0.99228 0.99070 0.00049 



The wavelengths included in the table are the observed wave- 
lengths in the frame of the planet (see text for additional de- 
tails). The time stamps included here are for the times at mid- 
exposure. The flux ratios are presented both before (raw) and 
after (corrected) EPD was applied. The flux ratios have also been 
normalized to the weighted mean OOT flux ratio (see Table [6] 
and text for more details). The full table is available online, and 
a portion is shown here so the reader can see the formatting of 
the table. 

and bluest off-line light curves are comparable, but the red 
off-line flux ratios (particularly for the reddest light curve) 
clearly lie at slightly higher values compared to the on-line 
flux ratios, indicating a smaller apparent planetary radius 
at those wavelengths. 

Ideally, when one has access to either a complete or 
partial transit light curve and baseline data acquired im- 
mediately before or after the transit event, one can fit a 
model to the data and estimate the transit depth from the 
model results. Due to the very long duration of HD 80606b's 
transit, we were not able to acquire baseline data on the 
night of the transit, thereby making this type of analysis 
impractical. However, thanks to several recent campaigns to 
observe a complete transit of HD 80606b and establish ac- 
curate orbital and p hysical parameters for this system via 
light curve modeling (|Winn et al.ll2009l ; iHebrard et alfcoid ; 
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Figure 8. Histograms of normalized flux ratios from the bot- 
tom of the transit light curve as shown in Figure [4}>. The his- 
tograms were generated using a bin size of 0.5 mmag. Each panel 
compares the on-line flux ratios with the off-line flux ratios. In 
each panel, the black (solid) histograms represent the 769.91-nm 
(on-line) light curve. The blue (dotted), brown (dashed) and red 
(dot-dashed) histograms are for the 768.76, 773.66 and 777.36 nm 
light curves and are shown in panels (a), (b) and (c), respectively. 
Panel (d) shows the histograms for all four wavelengths for further 
comparison. 

Iffidas et al1l20ld ; lshporer et alj|20ich . we do not need to fit 
a model to our partial light curve to achieve the goals of 
this paper. Instead, we consider only the middle ~4 hours 
of the transit light curve in our analysis (compared to the full 
duration of the bottom of the transit, which is ~6 hours), 
thereby minimizing systematic effects of stellar limb darken- 
ing (LD) as the strongest LD occurs during ingress, egress, 
and right after/before ingress/egress. Further, since we do 
not know the LD model for this star to the precision of our 
observations, adding such a model would not be useful for 
this study. Thus, we assume that LD is the same over all our 
bandpasses and that the transit ephemeris, impact parame- 
ter and transit duration do not vary with wavelength. The 
only parameter of which we assume changes with wavelength 
is the apparent planet radius {Hp). 

To investigate how the apparent planet radius changes 
with wavelength, we simply compute the weighted mean 
in-transit flux ratio [< SF/F >, which is proportional to 
the planet-star radius ratio, {R p /Ri,) 2 ] and its uncertainty 
for each wavelength. Specifically, we compute the weighted 
mean as 

n 

y^ y w,Fj 

< 5F/F >= i=l (1) 

where the weights, Wi, are equal to 1/(/3<t,:) 2 . Here, cy is 
the estimated photometric uncertainty weighted by some 



Figure 9. Standard deviation of the time-binned flux ratio mea- 
surements from the bottom of the transit (e.g. as shown in Figure 
[4f>) as a function of the number of data points per bin {N). Panels 
(a), (b), (c) and (d) show the standard deviations for the binned 
769.91, 768.76, 773.66 and 777.36 nm light curves. The amount 
of binning that could be performed varies for each light curve 
since the different wavelengths were observed a different number 
of times in a given observing sequence (see £|2,ll for details). The 
solid line in each panel represents the trend expected for pure 
white Gaussian noise (~ TV -1 / 2 ), normalized to the unbinned 
standard devation measured in our data. The dotted lines repre- 
sent the trend for Gaussian noise when normalized to the theoret- 
ical noise for our observations. The dashed curves are models fit 
to the standard deviation that include both white and red noise. 
The effect of red noise is obvious in all bandpasses. 

wavelength-specific factor (/}) in order to account for the 
presence of any red noise in each individual bandpass. 

To illustrate the effect of red noise on our measurements 
and the need for a re-weighting factor, the standard devia- 
tions (<tjv) of the in-transit and OOT time-binned flux ratios 
are shown in Figures [9] and [10] as a function of binning fac- 
tor {N) for each bandpass. The theoretical trend expected 
for white Gaussian noise (~ TV -1 / 2 ) is plotted as a solid 
curve, and we can see that for the in-transit data, the RMS 
deviates from the theoretical curve at large binning factors, 
indicating that red noise is present in most bandpasses (be- 
ing the least significant in the bluest bandpass). However, 
for the OOT data, our photometry appears to be generally 
consistent with the photon limit (although the bluest light 
curve suffers from small number stat istics) . 

Following me thods used by, e.g., IPont et all (|2006h and 
I Winn et all \200H ), we calculated explicit estimates for both 
the white {a w ) and red (ay) noise in each bandpass by solv- 
ing the following system of equations: 

<J\ = Ow + °"r (2) 
^2 

2 | 2 /o\ 

°~N = ~jy + °r- {&) 

The re-weighting factor, f3, is then computed as 
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Figure 10. Standard deviation of the time-binned OOT flux ra- 
tio measurements from April 2010 (e.g. as shown in Figure [6})) 
as a function of the number of data points per bin (N). Panels 
(a), (b), (c) and (d) show the standard deviations for the binned 
770.00, 768.86, 773.76 and 777.45 nm light curves. The solid line 
in each panel represents the trend expected for pure white Gaus- 
sian noise (~ iV" 1 / 2 ). The dotted lines represent the trend for 
Gaussian noise when normalized to the theoretical noise for our 
observations. The dashed curves are models fit to the standard 
deviation that include both white and red noise. Compared to the 
in-transit observations, red noise has a very minimal effect here. 
Deviations below the curve are likely due to small number statis- 
tics. These results demonstrate that narrow-band ground-based 
observations can provide very high-precision differential photom- 
etry. For a given bandpass, the combined precision ex ceeds that 
of Sp itzer llHebrard et al.|[201fj) or HST observations jPont et al.l 
2008). To the best of our knowledge, these represent the highest 
precision photometry for a 1.2-nm bandpass for ground or space 
observations. 



Figure 11. Normalized weighted mean in-transit flux ratio ver- 
sus observed wavelength (in the frame of the planet). The open 
triangles represent the flux ratios as computed for each light curve 
described in [[2] and [3] The solid circles represent the flux ratios 
computed after excluding outlying absolute flux values for each 
star from the analysis (see £|3 - IP ■ Note that the solid circles have 
been offset by 0.25-nm for clarity. The vertical error bars include 
a factor to account for the effects of red noise in both the in- 
and out-of-transit data. The "error bars" in the horizontal di- 
rection indicate the FWHM of each bandpass. The solid squares 
represent the mean in-transit flux ratios estimated from limb- 
darkened transit light curve models for HD 80606b. The lines 
show the predictions of planetary atmosphere models (see i|4,4l 
for more details). The inset figure shows the atmosphere models 
on a small vertical scale. While LD or night-to-night variability 
(of Earth's atmosphere or either star) could affect the overall nor- 
malization, the observed change in the flux ratio with wavelength 
is robust. 



Vr/icw/vN). Based on our fits to the red and white noise, 
we computed a re-weighting factor for each bandpass and 
applied it as stated above. We imposed a minimum value 
for P of 1, particularly for cases where red noise was negli- 
gible. 

The uncertainties for the OOT flux ratio are also 
weighted by the flux ratio, Fi, since two different exposure 
times were used during the OOT observations. Finally, the 
uncertainty on the weighted mean is computed as 



T <5F/F> 



1 



(4) 



We include the uncertainty on the weighted mean OOT flux 
ratio in our calculation of the mean normalized in-transit 
flux ratio and its uncertainty. The resulting spectrum of HD 
80606b (the normalized weighted mean in-transit flux ra- 
tios as a function of wavelength) is shown in Figure ITT1 and 
it clearly illustrates a difference between the flux ratios for 
the bluest bandpasses and those for the reddest bandpasses. 



While we find no significant difference between the flux ra- 
tios measured at 768.76 and 769.91 nm, we measure differ- 
ences of 3.02 ± 3.02 x 10" 4 and 8.09 ± 2.88 x 10~ 4 between 
observations at 769.91 and 773.66 and 777.36 nm. 

We list the weighted mean in-transit flux ratios (nor- 
malized by the weighted mean OOT flux ratios) as well as 
the weighted mean OOT flux ratios and their uncertainties 
in Table [6] In this table, we also include our fits to the white 
and red noise, as well as our estimates for (3. When calculat- 
ing the normalized in-transit flux ratio and its uncertainty, 
we also included the re-weighted uncertainty for the mean 
OOT flux ratio in our calculation. The error bars for the 
flux ratios given in Table [6] and shown in Figure [11] also take 
red noise into account. 



3.1 Effects of Earth's Atmosphere 

We consider the effect of random atmospheric variations 
(e.g., clouds) during the night of the transit as well as during 
the April baseline observations. As mentioned in gj2.2l large 
variations in the absolute flux of both the target and ref- 
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Table 6. Time- Averaged Flux Ratios and Noise Estimates 



X E (nm) 


\p (nm) 


A 5 (nm) 


< 8F/F > 


< J <SF/F> 












/3 


In- Transit 


768.60 


768.76 


768.60 


0.9901486 


1.32 xl0~ 4 


5.55 xlO" 


-4 


1.34 


xlO 


-4 


2.81 


769.75 


769.91 


769.75 


0.9901971 


2.08 xl0~ 4 


6.21 xlO" 


-4 


2.41 


xlO 


-4 


7.86 


773.50 


773.66 


773.50 


0.9904995 


2.19 xlO -4 


5.84 xlO" 


-4 


2.44 


xlO 


-4 


4.90 


777.20 


777.36 


777.20 


0.9910061 


1.99 xl0~ 4 


6.15 XlO" 


-4 


2.47 


xlO 


-4 


6.43 


Out-of- Transit 


768.86 




768.79 


1.1201531 


8.85 xlO" 5 


5.62 xlO" 


-4 


1.24 


XlO" 


-9 


1.00 


770.00 




769.93 


1.1200131 


5.05 xl0~ 5 


6.05 xlO" 


-4 


6.71 


xlO" 


■11 


1.00 


773.76 




773.69 


1.1202807 


8.46 xl0~ 5 


5.54 xlO" 


-4 


4.65 


xlO 


-9 


1.00 


777.45 




777.38 


1.1196202 


8.18 xl0~ 5 


5.78 xlO" 


-4 


8.29 


xlO 


-5 


1.45 



Xe is the observed wavelength from the GTC (i.e., from the Earth), Xp is the observed wavelength in the 
frame of the planet, and Ag is the observed wavelength in the frame of the star. Values for Xp are not 
given for the out-of-transit observations, as the planet was not transiting and was therefore not technically 
observed. The in-transit ratios refer to the relative flux ratio between the target and reference that has been 
normalized to the weighted mean OOT flux ratios (given at the bottom of the table). 



erence were observed towards the beginning and the end of 
the transit observations, with a few large fluctuations around 
the middle of the observations as well. Thus, to check if our 
measured in-transit flux ratios were affected by these fluctu- 
ations, we computed the weighted mean in-transit flux ra- 
tio for each bandpass after excluding outlying absolute flux 
measurements from our analysis. We specifically excluded 
any points that were greater than 3<r away from the mean 
of the flattest part of the spectrum measured for each band- 
pass and each star. After excluding outlying points from 
both the in-transit and April baseline data, we found that 
the new spectrum for HD 80606b shows a very similar shape 
as the original spectrum, albeit with the flux ratio in the red- 
dest bandpass differing the most from the original spectrum. 
However, we still measure a significant difference between 
the flux ratios in the on-line and reddest bandpasses. These 
results are included in Table [7] and shown in Figure [11] as 
the solid circles. 

3.2 Limb Darkening Effects 

So far our analysis has assumed that LD is the same between 
our different bandpasses, so LD should not affect the mean 
flux ratios for each bandpass differently. However, in prin- 
ciple, there is also the possibility that LD coefficients vary 
significantly in and out of narrow spectral lines. To inves- 
tigate the possibility that our spectrum's signature is a re- 
sult of our probing in and out of HD 80606's stellar spectral 
lines, we have computed quadratic LD coefficients for each of 
our ba ndpasses for a grid of stellar models [using PHOEBE; 
IPrsa fc Zwitterl (|2005f )]. We then generated theoretical limb- 
darkened light curves for each bandpass usin g the standard 
planet transit model of lMandel fc Agoll (|2002l ). We used stel- 
lar paramet e rs and uncertainties for HD 80606 as given by 
IWinn et all |2009) to estimate a range of LD coefficients 
to use in our models. We also input planet ary parameters 
and u ncertainties for HD 80606b as given bv lHebrard et all 
(2010). After computing light curve models for different 
combinations of LD coefficients and planetary parameters, 
we computed the mean model flux ratio over the bottom 



of each transit light curve (the 4 hours centered around 
mid-transit). We include the resulting model spectrum in 
Figure [TT] as solid squares. This particular spectrum was 
computed based on using a median set of LD coefficients, 
but all the model results were similar over the range of LD 
coefficients used. The median linear and quadratic LD coeffi- 
cients (wi,u 2 ) are (0.392,0.229), (0.388,0.233), (0.391,0.230) 
and (0.376,0236) for the 768.76, 769.91, 773.66 and 777.36 
nm bandpasses. 

While small differences in LD exist between the differ- 
ent bandpasses, the mean model flux ratios differed by only 
a very small amount (< 2 x 10 -5 ) between the different 
bandpasses. From this, we conclude that LD is most likely 
not the cause of the large variations in our observed spec- 
trum. However, we note that PHOEBE (as well as other LD 
codes) has not been calibrated in and out of narrow spectral 
lines. We also note that the models show that the bottom 
of the light curve is in fact not flat due to LD. However, 
based on our calculation of the mean model flux ratio over 
the limb-darkened transit bottom for each bandpass, this 
should not affect the magnitude of the variations we mea- 
sure in our observed spectrum. Due to LD effects, the overall 
normalization of the spectrum may be affected. 

3.3 Transit Colour 

In Figure [12] we present the colour of the normalized in- 
transit flux ratios, computed by dividing each point in the 
off-line bandpasses by the average of each pair of on-line 
points around those off-line points. We find that the colour 
between the bluest bandpass and the on-line bandpass is 
consistent with zero, with a mean value of 6.30L6.04 x 10~ J 
(computed following the method described in fj3]). The 
mean colour of the 773.66 nm and on-line bandpasses is 
-3.57L0.63 x 10" 4 , and the mean colour between the red- 
dest and on-line bandpasses is — 8.99L0.62 x 10 -4 . 

We also present the standard deviation of each colour 
for a number of binning factors in Figure [T3] We find that 
the trend for each colour is consistent with having only white 
noise in each of our colours. This is also confirmed by fitting 
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Table 7. Time- Averaged Flux Ratios and Noise Estimates (Outlying Absolute Fluxes Excluded) 



X E (nm) Ap (nm) X s (nm) <8F/F> o" <<SF / F> a w a r j3 



In- Transit 



768.60 


768.76 


768.60 


0.9901021 


4.83 


XlO" 


5 


4.88 


xlO" 


-4 


3.90 


xitr 


5 


1.00 


769.75 


769.91 


769.75 


0.9901218 


1.76 


XlO" 


4 


6.29 


xHT 


-4 


2.03 


xlO" 


4 


6.44 


773.50 


773.66 


773.50 


0.9904292 


2.15 


xlO" 


-4 


5.82 


xlO" 


-4 


2.40 


xlO" 


4 


4.80 


777.20 


777.36 


777.20 


0.9907548 


2.42 


xlO" 


-4 


5.89 


xlO" 


-4 


2.90 


xlO" 


4 


7.88 



Out-of- Transit 



768.86 




768.79 


1.1202527 


1.01 


xlO" 


4 


5.66 xlO" 


-4 


7.78 


xlO" 


9 


1.00 


770.00 




769.93 


1.1201327 


1.20 


XlO" 


4 


5.89 XlO" 


-4 


1.13 


XlO- 


4 


2.19 


773.76 




773.69 


1.1203722 


8.98 


Xl0" 


5 


5.77 XlO" 


-4 


1.35 


xlO" 


-9 


1.00 


777.45 




777.38 


1.1196849 


5.99 


Xl0" 


5 


5.82 xlO" 


-4 


3.36 


xlO" 


-9 


1.00 



Same as in Table [6] but the flux ratios listed here are those computed after excluding outlying absolute flux 
measurements from the analysis. 




Figure 12. Colours of the normalized in-transit flux ratios. The 
different panels show the colour as computed between each off- 
line bandpass and the on-line bandpass (after binning the on-line 
data to the number of points in each of the off-line bandpasses). 
The dashed line in each panel illustrates where the colour equals 
zero. The data has not been explicitly offset, and that there are 
no obvious systematics seen in any of the colours. 



the white and red noise explicitly for each colour. Consider- 
ing that the red noise is estimated to be less than ~lxl0 -8 
for each colour, white noise clearly dominates the uncertain- 
ties in the transit colour. 

As explored in £13.11 we also compute mean colours af- 



Figure 13. Standard deviation of the time-binned colour mea- 
surements from the bottom of the transit (as shown in the dif- 
ferent panels in Figure 1126 . The different panels show the stan- 
dard deviations for the different colours as presented in the pan- 
els in Figure IT2l with panels (a), (b) and (c) respectively show- 
ing the standard deviations for the 768.76 - 769.91 nm, 773.66 
- 769.91 nm and 777.36 - 769.91 nm colours. The solid line in 
each panel represents the trend expected for pure white Gaussian 
noise (~ N~ V 2 ). The dotted lines represent the trend expected 
for Gaussian noise when normalized to the unbinned theoretical 
uncertainties for these observations. There is no obvious presence 
of red noise at large binning factors. 
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0.0000 




770 772 774 776 

Off — Line Wavelength (nm) 

Figure 14. Mean colour of the in-transit flux ratios as com- 
puted between each off-line bandpass and the on-line bandpass. 
The open triangles represent the colours as computed in i|3,3l and 
illustrated in Figure 1121 The solid circles represent the colours 
computed after excluding outlying absolute flux measurements 
for each star from the analysis (see i|3.1l l. The errors bars repre- 
sent the lo" uncertainties. The dashed line illustrates where the 
colour equals zero. We arbitrarily set this point equivalent to an 
apparent planet radius of 1 (i.e., we let the measured radius in the 
on-line bandpass be the baseline radius of HD 80606b). The mean 
colours around the 773.66-nm bandpass are essentially equal for 
both sets of points, so the two data points appear as one. 



atmosphere is real and that the differences in the flux ratio 
are significant. 



4 DISCUSSION 

4.1 Interpretation of Light Curve Shape 

First, we compare our light curve (integrated over all 
bandpasses) to sim ultaneous observations from Spitzer 



Hcbrard et al.l 



2010) and other grou nd-based observatories 



Shporer et alj|20ich . In particular, iHebrard et~ail (|2010h 
identified a bump in the in-transit light curve that oc- 
curred within the hour before their estimated time of mid- 
transit and pondered whether it could be due to a Moon 
or spot crossing. Under the Moon hypothesis, the magni- 
tude of the bump should be wavelength independent. If the 
bump were due to a spot, then one would expect an even 
greater feature in the optical. We do not find any evidence 
for a coincid ent bump [regardless of whether we adopt th e 
ephemeris of IHebrard et all (|201Ch or lShporer et al.l l|2010l )]. 
Thus, the bump is unlikely to be due to either a Moon or 
starspot. If anything, we find possible evidence of a bump 
occurring after mid-transit, b ut this feature was not de- 
tected by IHebrard et all (|2010l ). If we assume our candidate 
bump is not a result of instrumental systematics, and we 
compare our candidate bump as observed in the different 
wavebands, we find that the size of the putative bump is 
smallest in the bluest bandpass, providing further evidence 
against a starspot. Furthermore, since the magnitude of the 
bump varies slightly for each bandpass, this provides addi- 
tional evidence against the existence of a Moon. Future high- 
precision, multi- wavelength observations could help provide 
additional constraints on the light curve shape. 



ter excluding outlying absolute flux measurements from our 
analysis. After excluding those data points, we estimate the 
mean colours between each off-line and the on-line band- 
passes to be 1.79±6.60 x 10~ 5 , -3.54±0.62 x 10~ 4 , and 
-6.92±0.54 x 10" 4 (from bluest to reddest). Both these 
mean colours and those discussed above are plotted in Figure 
1141 The colours are comparable between the two data sets, 
with the colour of the reddest bandpass having the only mea- 
surable difference between the two sets. Furthermore, Figure 
1141 illustrates that not only is there a significant change in 
the colour during transit, but also that the magnitude of 
the change is equivalent to a large change in the apparent 
planet radius. At the reddest wavelengths, we clearly mea- 
sure a change of over 3% (and as much as 4.2%, based on the 
flux ratios that do not exclude outlying absolute flux mea- 
surements) in the apparent radius of the planet compared 
to the planet's apparent radius in the on-line bandpass. 

Overall, these colours agree with the magnitude and di- 
rection of the differences measured between the weighted 
mean in-transit flux ratios for the different bandpasses (see 
Furthermore, the differences between the colour of 
the bluest-to-on-line bandpasses and the reddest-to-on-line 
bandpasses has greater than 5er significance. Since our colour 
measurements match the magnitude and direction of the dif- 
ferences in the flux ratios as measured from our spectrum, 
we conclude that our measured spectrum of HD 80606b's 



4.2 Comparison to Previous Observations 

Next, we note that our measured in-t ransit flux ratios diffe r 
slightly from the flux ratio given by IHebrard et ail l|2010h . 
This is at least partly due to the different bandpasses used. 
There could also be a systematic uncertainty in the overall 
normalization of our transit depths. If our goal had been 
to measure the transit depth precisely, we would have re- 
quired observations taken just before and after the transit 
event. In this case, ground- based observations spanning the 
full transit were not feasible due to the extremely long tran- 
sit duration. Thus, we normalized our in-transit light curves 
by OOT observations taken on a different night. While our 
observations resulted in a very high precision for differen- 
tial measurements of the transit depth in each bandpass, a 
change in the observing conditions between nights could re- 
sult in the transit depths all being affected by a common 
scaling factor. 

To confirm that the change in the apparent planetary 
radius with wavelength is based on a robust estimate of the 
OOT flux ratio despite using baseline observations separated 
by four months from the transit observations, we estimated 
the weighted mean in-transit flux ratios as before, but nor- 
malized them against the lower quality OOT data taken on 
2010 January 15. For reference, we include the results of the 
aperture photometry for this data set in Table [8] and the 
flux ratios before and after EPD in Table [9] We found that 
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Table 8. Absolute OOT Photometry from 2010 January 15 

A (nm) HJD F targe t F ref 

768.60 2455211.7489 1077950 962913 

769.75 2455211.7487 1210710 1082178 



773.50 2455211.7492 1020447 911377 



777.20 2455211.7495 1301394 1161982 



Columns are similar to Table [T] except the wavelengths included 
in the table are the wavelengths as observed from the GTC (see 
text for additional details). The full table is available online, and 
a portion is shown here so the reader can see the formatting of 
the table. 

Table 9. Relative OOT Photometry from 2010 January 15 



A (nm) HJD F ratio F ratio Uncertainty 

(raw) (corrected) 

768.60 2455211.7489 1.11947 1.12007 0.00128 



769.75 2455211.7487 1.11877 1.11775 0.00121 



773.50 2455211.7492 1.11968 1.11940 0.00131 



777.20 2455211.7495 1.11998 1.11910 0.00117 



4.3 Lack of a K I Line Core 

As illustrated in Figure [TT1 there is no significant difference 
between the observations acquired in the core of the K I 
line and slightly to the blue. Given the 1.2-nm FWHM, a 
Doppler shift of > 200 km s _1 would be needed to shift the 
line core out of the on-line bandpass. This is greater than 
the escape speed from HD 80606b (~121 km s" 1 ). Thus, 
we place a 3a limit on the strength of the K I line core of 
3 x 10~ 4 (for our 1.2-nm FWHM bandpass). 

By itself, the lack of a line core is most naturally ex- 
plained by a lack of K I at the altitudes probed by trans- 
mission spectrophotometry. This could occur if: 1) there is a 
significant bulk under-abundance of potassium, 2) the potas- 
sium has condensed into clouds and/or molecules, 3) there 
is a cloud or haze layer above the region capable of causing 
significant potassium absorption, and/or 4 ) the potassium 
has been photoionized l|Fortnev et alJlibost ). In the previous 
case of HD 209568b, theoretical investigations of the unex- 
pectedly weak Na I absorption showed that the observed 
feature dep th is particularly se nsitive to the extent of cloud 
formation dFortnev et all 120031 ). In the case of HD 80606b, 
the highly eccentric orbit results in flash heating near peri- 
center and extreme temperature variations over the orbital 
period. At the time of transit, the star-planet separation 
is ~0.3 AU, so the equilibrium temperature is ~ 500 K. 
Based on Spitzer observations, cooling is sufficiently rapid 
that the planet i s expected to have co oled between percen- 
ter and transit IjLaughlin et al.l 120091 ). Thus, both sodium 
and potassium are predicted to have condensed into clouds. 
Thus, we conclude that the lack of a K I core could easily be 
due to potassium having condensed into clouds before the 
time of transit. 



The columns are similar to Table \E\ The full table is available 
online, and a portion is shown here so the reader can see the 
formatting of the table. 

despite the large scatter in that OOT data, the normalized 
in-transit flux ratio (and therefore, the apparent radius of 
the planet) still changes significantly with wavelength and 
maintains the same shape as shown in Figure 1111 We con- 
clude that the large change in transit depth from the 768.76 
and 769.91 nm bandpasses to the 773.66 and 777.36 nm 
bandpasses is a robust result. We also emphasize that the 
atmosphere was much more stable during the April observa- 
tions than the January baseline observations, so we still rely 
on the April baseline observations for our primary analysis. 

We also tested whether our results were sensitive to the 
aperture radius used for photometry. We tried a variety of 
annuli for the apertures for both the in and out-of-transit 
data sets. In all cases, we see the trend of increasing flux ra- 
tio with wavelength and found very similar results to those 
presented here. The only difference occurred for the largest 
apertures, in which case the weighted mean in-transit flux 
ratio on the K I feature is slightly smaller than the flux 
ratio at the bluest wavelength. Even for this choice of aper- 
tures, the fluxes in the 773.66 and 777.36 nm bandpasses are 
not significantly different (though the error bars are slightly 
larger) . 



4.4 Planetary Atmosphere Models 

In an attempt to model our observations, we c onsidered both 
a con ventional 1-D "cold" atmosphere model (|Fortnev et al.l 
boicf ) (solid line in Fig. Ill[) and a similar model, but with ar- 
bitrary additional heating to raise the effective temperature 
by 500 K (dotted line in Fig. Ill|) . Both mo dels have been 
norma lized to the stellar radius estimated bv lHebrard et all 
(2010), and assume a star-planet separation of 0.3 AU (i.e., 
the distance between the star and planet when the planet 
transits). Chemical equilibrium and a standard pressure- 
temperature profile for HD 80606b are assumed. In the 
"cold" atmosphere model, the planet's (apparent) radius 
at 10 bar was adjust ed to match the radius measured by 
Hebrard et al l (|2010l ) at 4.5 microns. In the "hot" atmo- 
sphere model, the temperatures in the upper atmosphere 
range from ~300— 500 K, even with the additional heating. 
The higher temperature increases the observed planetary ra- 
dius at all wavelengths, and slightly increases the peak to 
trough distance of the features, but t he planet's rad i us wa s 
not adjusted to match the radius from lHebrard et alj (|2010l ). 
At these temperatures, most of the potassium is expected to 
have formed condensates, significantly reducing the K I ab- 
sorption feature. As the inset in Figure[TT]illustrates, neither 
the "hot" and "cold" atmosphere models predict a signifi- 
cant feature due to K I absorption. 
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4.5 Change in Apparent Radius with Wavelength 

While we do not detect the K I core, we find relatively large 
differences (3.57 ± 0.63 x 1(T 4 and 8.99 ± 0.62 x 10 -4 ) be- 
tween the colours of the on-line bandpass and the band- 
passes to the red (773.66 and 777.36 nm). Clouds and hazes 
would suppress both the core and wings of the absorption 
feature. A similar observation for a typical hot- Jupiter could 
be readily interpreted as strong absorption in the wings of 
the potassium line due to absorption by pressure broadened 
potassium at lower altitudes, while potassium at h igher al- 
titudes has been photoionized (jFortnev et alj|2003f ). 

However, in our observations, the magnitude of the dif- 
ference in absorption at the two blue and two red wave- 
lengths appears too large for such a model. One could ex- 
pect such observations to probe the lower atmosphere over 
~10 scale heights (H), from a pressure of ~100 mbar to 
~1 microbar. Assuming the planet has reached a thermal 
equilibrium for the star-planet distance at the time of tran- 
sit and a 500K upper atmosphere temperature, the scale 
height would be H ~ 20 km. Thus, one might expect to see 
changes in the apparent radius of the planet on the order of 
~ 200 km. Our observations suggest a much larger change 
in the apparent radii (up to ~ 4.2% or ~ 2900 km) when 
comparing observations in the K I line core and the reddest 
bandpass. The scenario described above would suggest that 
these observations probed ~145 scale heights in the atmo- 
sphere of HD 80606b, or pressures of less than ~ 10 -55 bars, 
which is well into the exosphere. Such a large number of scale 
heights is not realistic, implying that the absorption is orig- 
inating from a part of the atmosphere much hotter than 500 
K. Fortunately, the temperature is expected to rise r apidly 
to th ousands of Kelvin above one planetary radius |Yellel 
l2004h . 



exosphere does not exist. We hope that our observations 
will stimulate theoretical models for the observable effects 
of exoplanet exospheres on transmission spectroscopy and 
spectrophotometry. 

Given the planet's high surface gravity and any reason- 
able choice of planetary parameters, a ~ 4.2% change in the 
planet's apparent radius requires a very dramatic change in 
the pressure at which the slant optical depth reaches unity, 
between 770 and 777 nm. Thus, we conclude that absorption 
at high altitude and temperature is the most likely explana- 
tion for the large change in the apparent planet radius. 

4.7 Possibility of Other Absorbers 

Next, we consider whether another absorber might be re- 
sponsible for the observed change in apparent planet ra- 
dius. Methane can be active in this region of the spectrum. 
However, methane would be unstable at the high tempera- 
tures of an exosphere or wind. Both of the models in £|4 .41 
include methane at all temperatures at which it would be 
stable, around <1000 K. In the wavelength regime that we 
observed, the opacity of methane is largest at 778 nm and 
smallest at 769 nm, so its presence would produce the op- 
posite trend from what is shown in the data. 

The observed wavelengths were also chosen to avoid wa- 
ter vapor (which is also unstable at high temperatures). We 
are not of aware of any other absorber which could explain 
the large change in apparent planet radius, and consider K I 
the most likely absorber. Nevertheless, we cannot rule out 
the possibility that HD 80606b 's exosphere possesses an ab- 
sorber that is something other than K I on account of the 
incomplete opacity database. 



4.6 Absorption by an Exosphere 



Based on the model of £14.41 we would estimate that our ob- 
servations have probed ~145 scale heights in the atmosphere 
of HD 80606b, or a pressure of less than 10 -55 bars. How- 
ever, t hese estimat es assume an atmospheric temperature of 
500 K. lYellel i|2004l ) finds a steep rise in the temperature from 
~350 to 10000 K from 1 R p to 1.1 R p for a planet at 0.1 AU 
from the Sun. If we use their model as a rough guideline, and 
if we assume a temperature of 2000 K between 1 and ~1.04 
R p for HD 80606b, the 2900 km measured change in the 
apparent radius would imply that the observations probed 
~36 scale heights, or to a pressure of less than 10 -14 bars. 
Regardless of whether we assume a temperature of 500 K 
or 2000 K, the implied pressures are indicative of those that 
would exist in an exosphere. 

The models and opacity database of £14.41 are not 
complete for the temperature and pressures of the exo- 
sphere. The opacity database used extends to temperatures 
of ~2600 K and ~1 microbar and is not intended to 
describe opacity so u rces in an e x osphere or win d (e.g. , 
Vidal-Madiar et all 1200 j . 120081: IrBailester et al.l 120071; 
Ehrenreich et all 120081; iLecavelier Pes Etangs et al.l 



200; 



^ecavelier des Etangs et all 
Ben-Jaffel fc Sona Hosseinil l201ol l To the best of our 
knowledge, an exospheric model that predicts the loca- 
tion and strength of absorption features arising from the 



4.8 Absorption by a Wind 

If the ~ 4.2% change in the apparent radius is due to absorp- 
tion by K I at high altitude, then it is not obvious why the 
observations on the K I core (769.91 nm) are not significantly 
different from the observations slightly to the blue (768.76 
nm). One possibility is that the line core was shifted out 
of the on-line bandpass. Given the ~ 1.2-nm FWHM, this 
would require a Doppler shift of > 200 km s" 1 . A blue shift 
of 225-km would place the core halfway between the 768.76 
and 769.91 nm bandpasses. A somewhat smaller Doppler 
shift plus Doppler broadening might also reduce the sig- 
nal strength. In any case, the velocities required would be 
greater than the escape speed from HD 80606b (~121 km 

While a velocity exceeding the escape speed is some- 
what concerning, it is not out of the question for a wind 
being driven from the exosphere. In fact, similar observa- 
tions of other planets also appear to find an unexpectedly 
large Doppler shift. Speci fically, a large b l ue shi ft has been 
found in all cases. E.g., iRedfield et alj (|2008T ) found an 
unexpected blue shift of the c o re of the Na I absorption 
for HD 189733b. ISnellen et all l|2010ri detected a 2 km s _1 
blueshift in the upper atmosphe re of HD 209458 b with 
observations of CO. Additionally, iHolmstrom et all (|2008T ) 
reported Lyman-a absorption around HD 209458b at 
wavelength offsets corresponding to velocities of several 100 
km s^ 1 , but there was no information about the Lyman-a 
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core as it is not observable due to Earth's geo-corona. Much 
like our observations of HD 80606b, there is considerable 
uncertainty regarding the ori gin of the absorption and 
Dopp ler shift f or HD 209458b |Lecavelier Pes Etangs et al.l 



200; 



Lccavclicr dcs Etangs et al. 



20101 ; 



Ben-Jaffel fc Sona Hosseinil l2010h . Proposed mecha- 



nisms include radi ation pressure and interaction with a 
stellar wind (e.g., iTian et ail 120051; iGarci'a Mufio3 120071 ; 
iMurrav-Clav et aLll2009l ; lEkenback et alJl201Cf ). and in par- 
ticular we note that models of HD 209458b's atmosphere 
match observations better if it is assumed that the line 
core is obscured. Our observations could be explained if 
a similar mechanism operates on HD 80606b and heavy 
elements (i.e., potassium) are mixed into the wind. 

In the case of HD 80606b, the dynamics of the exosphere 
and any planetary wind is almost certainly quite complex. 
The planet has the largest semi-major axis of any confirmed 
transiting planet (0.455 AU), but it follows such a highly 
eccentric orbit (e = 0.93) that the star-planet separation of 
HD 80606b at periastron is ~2/3 that of HD 209458b. Thus, 
HD 80606b experiences strong and rapid heating of the at- 
mosphere near pericenter. The large and rapid changes in 
the incident stellar flux and temperature, as well as the stel- 
lar wind flux could lead to episodic mass loss following each 
pericenter passage dLaughlin et al.ll2009h . Based on the ob- 
served X-ray flux (iKashvap et al, I l2008h and mass loss corre- 
lation l|Wood et al.ll2005l ). HD 80606 could have a mass loss 
rate as much as ~100x stronger than HD 209458, provid- 
ing a much stronger stellar wind to drive a wind from HD 
80606b. The rapid contraction and expansion of the Roche 
lobe around each pericenter could further complicate the 
dynamics of the exosphere and planetary wind. 



4.9 Potential Systematics 

4-9.1 Excluding Telluric Absorption 

The usual suspect in ground-based observations is variabil- 
ity in the telluric absorption. At our observed wavelengths 
there is very little absorption. The only two species that con- 
tribute any appreciable absorption are water and oxygen. In 
particular, there is a lack of absorption from carbon dioxide 
or methane in our observed bandpasses. Oxygen is generally 
well mixed in the atmosphere. Thus, we expect any variabil- 
ity due to oxygen has been removed in our data reduction 
procedure, which normalizes each observation of HD 80606 
by the flux of HD 80607 taken at the same time and using 
the same bandpass. Thus, we rule oxygen absorption out as 
a potential systematic. 

Since water can be very anisotropically distributed in 
the atmosphere, one could worry that the 20 arcsec separa- 
tion between HD 80606 and HD 80607 might allow for vari- 
ations in the water absorption that are not removed by cali- 
bration. However, the two bandpasses to the red of K I were 
specifically chosen to be at wavelengths that avoid water 
absorption. Thus, even in the scenario that the on-line and 
blue bandpasses were contaminated by water absorption, we 
still measure a ~ 2.7% change in the apparent planet radius 
between the two reddest bandpasses (both of which should 
be substantially free of telluric absorption). From this, we 
conclude that our primary result of measuring a large change 



in the apparent radius with wavelength is not the result of 
variable telluric water absorption. 

However, in an effort to confidently rule variable tel- 
luric absorption out as a potential source for systematics, 
we construct an alternative model for the spectrum based on 
changing the level of water vapor absorption. Specifically, we 
integrated our bandpasses over high-resolution model trans- 
mission spectra for telluric water vapor and oxygen. Using 
the TERRASPEC code (Bender et al. 2011, in prep), we 
computed model transmission spectra for two different air- 
masses (representing the mean airmass over the transit bot- 
tom as observed in January and the mean airmass during 
the baseline data observed in April) and three different water 
vapor levels (1,5 and 10 mm). We then integrated our band- 
passes over each spectra and computed the relative trans- 
mission for the different bandpasses for every possible com- 
bination of water vapor towards HD 80606 and HD 80607. 
We integrated over the appropriate bandpasses for each set 
of observations, as the bandpasses used in April were cen- 
tered at slightly different wavelengths than for the January 
observations. Our goal was to determine if the transmission 
spectrum would have a similar signature as our observed 
spectrum if there was a difference in the water vapor to- 
wards HD 80606 and HD 80607 during either or both of the 
January and April observations. For example, we took the 
integrated transmission for a water vapor level of 10 mm 
(towards HD 80606) divided by the integrated transmission 
for a water vapor level of 1 mm (towards HD 80607) based 
on the mean airmass in January. Then, we divided that re- 
sult by a similar ratio based on the spectra for the April 
observations. We computed this ratio for all combinations 
of water vapor and compared the results. Realistically, the 
water vapor was most likely below 6 mm for both the Jan- 
uary a nd April observations [based on lGarcfa- Lorenzo et al.l 
(2010)], but we approach this issue with much caution and 
therefore discuss the results for the 10 mm water vapor level 
as well. 

From our results, we can make several arguments 
against variable water vapor absorption and/or the different 
wavelengths observed in January and April being the cause 
of our spectrum's signature. First, since our measurements 
are multiply differential (comparing the target to the refer- 
ence in-transit to the target to the reference out-of-transit), 
we minimize any such effects from our measurements. Sec- 
ond, even if the water vapor column towards HD 80606 and 
HD 80607 differed by an average of 10 mm on one of the 
nights, it would result in a difference of only 0.0058% (if 
the January water vapor differed by 10 mm) or 0.034% (if 
the April water vapor differed by 10 mm) between the red- 
dest and on-line bandpasses. This difference in transmission 
based on the wavelengths observed in April is less than half 
of the actual measured difference between the flux ratios in 
these two bandpasses. Further, the separation between HD 
80606 and HD 80607 is only 20 arcsec, so it is extremely un- 
likely that the time-averaged water column towards the two 
stars would differ by 10 mm. Further, an untenably large 
water column, inconsistent with the observational condi- 
tions, would be required to explain the observed difference 
of ~8xl0 -4 between the on-line and reddest bandpasses. 
Third, even if the average water column towards the two 
stars did differ by that much on one of the nights, the re- 
sulting colours differ from what we observe. I.e., the hy- 
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pothesis that our measurements are primarily due to atmo- 
spheric variability would predict the two bluest bandpasses 
to be comparable in some cases and differ largely in oth- 
ers, while the two reddest bandpasses are comparable in all 
cases. While we do observe the flux ratios in the two bluest 
bandpasses to be comparable, we see a significant difference 
between the flux ratios for the two reddest bandpasses. Fur- 
ther, the magnitude of the differences between the bluest 
and reddest bandpasses is observed to be much larger than 
what the difference would be if they were caused by variable 
atmospheric absorption. 

Thus, we estimate that for the four different band- 
passes, the effect of variable atmospheric absorption would 
be less than (0.0024%, 0.0015%, 0.0075%, 0.0073%) x [ 
< mm of H 2 towards HD 80606 > - < mm of H 2 to- 
wards HD 80607 > ]/[10 mm of H 2 0] at the wavelengths 
and airmass observed at in January, or less than (0.004%, 
0.027%, 0.0041%, 0.0068%) multiplied by the same ratio 
given above at the wavelengths and airmass observed at in 
April. However, assuming that the difference in transmis- 
sion is neglible between the target and reference for both 
the January and April observations, then the fact that the 
stars were observed at different wavelengths and airmasses 
on those nights should be irrelevant. 

Finally, we note that spectrophotometry using a 
narrow-band tunable filter is much less prone to systematics 
than spectroscopic observations. The lack of a slit, the simul- 
taneous use of a very good reference star, rapid switching 
between bandpasses, and the multiply differential nature of 
our measurement should all minimize the effects of telluric 
absorption. While the OH lines are variable, the sky sub- 
traction in our data reduction process removes the emission 
to a high degree of precision. Finally, we see no evidence, 
in our atmospheric transmission models, of absorbers that 
could account for the signal detected. 

To first order, the effects of atmospheric extinction are 
corrected by measuring flux ratios relative to HD 80607. We 
expect negligible second order differential extinction, since 
the target and reference stars are of the same spectral type 
and separated by only 20 arcsec. Since the magnitude of this 
effect scales as the square of the filter bandpass, our use of 
such a narrow bandpass further minimizes second order dif- 
ferential extinction, allowing this technique to be applied to 
other targets with reference stars that differ in temperature. 

We do not consider differential extinction to be a viable 
explanation for the effect seen in Figure [TT1 Nevertheless, we 
performed an additional check, in which we do not perform 
relative photometry between the target and reference. We 
compare the ratio of the absolute flux of the reference star 
in the reddest bandpass and the bandpass centered on K I 
as measured on the night of the transit to the same ratio 
as measured on the night the OOT observations were taken 
(2010 April 4). We estimate a ratio of 0.98498 ± 0.00093, 
equivalent to a colour deviation of ~1.5% between the two 
nights. This provides an upper bound on the effects of atmo- 
spheric variability, including differential extinction. The ac- 
curacy of our primary analysis should be considerably higher 
thanks to the use of relative photometry to correct for at- 
mospheric variability. 



4-9.2 Excluding Instrumental Effects 

With tunable filter imaging, the photons for each observed 
wavelength lands on the same pixels, eliminating concerns 
about spatial variations in the flat fielding. However, the 
normalization of the flux measurements is affected by the 
wavelength-dependence of the pixel sensitivity. To minimize 
this effect, we took dome flat fields for each observed wave- 
length and corrected the science frames taken at each wave- 
length with their respective flat fields. 

Furthermore, we guard against possible non-uniformity 
in the shutter motion, which could result in the system- 
atic effect of producing slightly different exposure times for 
the target and reference star, depending on where they are 
located on the CCD chip. This systematic effect is more no- 
ticeable for shorter exposure times, so we guard against it 
by following an observing sequence that repeats after seven 
exposures, so that the subsequent set of exposures occurs 
with the shutter motion in the opposite direction. 

Depending on the orientation of the tunable filter, the 
observed wavelength can drift due to the rotation of the in- 
strument during the observations. For our observations, the 
tunable filter was tuned before observations, in the middle 
of the transit and at the end of the observations. No drifts 
larger than 0.1-nm occurred. 

Finally, we conducted a thorough investigation into the 
possibility of saturation and/or non-linearity as a source of 
systematic effects. A majority of the peak counts during 
our observations were well below the saturation threshold 
(~65,000 ADUs), and for standard observing modes lin- 
earity is guaranteed up to ~65,000 ADUs, so non-linearity 
should not be an issue. However, as we use a non-standard 
observing mode on OSIRIS in order to read out the CCDs at 
the highest rate possible (and thereby greatly reduce dead 
time), it is worthwhile to investigate whether non-linearity is 
an issue. Therefore, we discuss here several checks for non- 
linearity, where we arbitrarily assume that 45,000 counts 
(~30,821 ADUs, based on the gain of 1.46 e" per ADU) is 
the level at which non-linearity might begin. 

First, we checked if the average number of counts from 
the flat-fields taken for each bandpass had a linear depen- 
dence with exposure time, as we had taken flat-fields at sev- 
eral different exposure times. We fit a line to all measure- 
ments of the mean flat-field counts (for 5 different exposure 
times), and then we fit another line to the data but excluded 
measurements that were near or above 45,000 counts. To see 
if including measurements at higher counts resulted in non- 
linearity, we compared the slopes and y-intercepts of the two 
best-fit lines. After comparing the best-fit solutions between 
the different bandpasses and for the different series of flats 
taken in January and April, we find that it is not obvious 
that any one set of flats displays significant non-linearity 
compared to the others. 

Second, we investigated fitting a quadratic function to 
the flat-field counts for both the January and April flat- 
fields. After comparing the best-fit coefficients for the dif- 
ferent bandpasses, we found that at least one set of coef- 
ficients deviated significantly from the coefficients for the 
other bandpasses. While there might be an obvious outlier 
in terms of one bandpass that might be affected by non- 
linearity for each set of flats, the supposed outlier is different 
for the January and April flats. We again conclude that it 
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is not obvious which, if any, of our bandpasses is displaying 
significant non-linearity. 

Third, we investigated the possibility of non-linearity by 
computing the colour (between each off-line bandpass and 
the on-line bandpass) and seeing how it varied with respect 
to the average on-line flux per pixel (estimated by dividing 
the total absolute on-line target flux by the target's FWHM 
squared). We computed the median colour for exposures 
where the average flux per pixel was below 45,000 and for 
exposures where the average counts were above 45,000. We 
then estimated the difference in the median colour for those 
two sets of exposures. We found that in the near-red band- 
pass (773.66-nm), non-linearity most likely does not play 
a role, as the median colours below and above the 45,000 
count level differ by an insignificant amount. However, in 
the bluest and reddest colours, we do see a slight correlation, 
with the median colours differing by comparable amounts. 
This is not what we might expect to see if non-linearity were 
causing systematic effects in our observed spectrum, since 
we do not observe a comparable colour difference in the in- 
transit flux ratio between the blue-on-line and reddest-on- 
line colours. 

We also computed the colour deviations for the April 
baseline data, and we found that the colours below and 
above the 45,000 count mark are slightly larger than the 
in-transit colour deviations (but these were computed us- 
ing a combined data set for two different exposure times, 
which could affect these estimates). Regardless, we still find 
that the smallest difference in the colours is in the near-red 
bandpass, and the differences are comparable for the bluest 
and reddest bandpasses, even though in both the observed 
in-transit and out-of-transit flux ratios we see the smallest 
difference in the flux ratio between the bluest bandpasses 
and the largest between the on-line and reddest bandpass. 

In summary, we conducted several checks for non- 
linearity. We conclude that any effects of non-linearity are 
either too insignificant to affect our photometry or they are 
not correlated with the data. 

4-9.3 Possible Non-Planetary Astrophysical Effects 

Lastly, we consider potential astrophysical systematics such 
as stellar variability. Observations in all four bandpasses 
were obtained during the same transit. If observations using 
different bandpasses had been made during different tran- 
sits, then the interpretation would be ambiguous, as stellar 
variability (e.g., spots that the planet does not necessarily 
pass over) could result in apparent changes in the in-transit 
flux ratio. For the large change in apparent radius to be due 
to stellar variability, there would need to be a ~ 4.2% change 
in the colour of either the target or reference star. Such large 
variability over a small r ange of wavelengths is a priori un- 
likely for solar-like stars (jHebrard et al.l 120101 ). However, as 
suggested by the referee, we have estimated how spotted HD 
80606 would have to be to cause a difference of ~8x 10 -4 in 
the flux ratios in the on-line and reddest bandpass. 

We computed the blackbody flux for HD 80606 (T eff 
~ 5572 K) then integrated the flux over each bandpass to 
estimate the total flux observed in each bandpass. We then 
completed similar calculations for a spot assuming a tem- 
perature 1000 K cooler than HD 80606 and a spot radius 
equal to the planet's radius. After computing the ratio of 



the integrated spot flux to the integrated star flux for some 
N spots, we found that about 26 spots with the above prop- 
erties would have to exist on the surface of HD 80606 during 
the transit observations in order to produce a difference in 
the on-line and reddest flux ratios of about 8xl0~ 4 . That 
is equivalent to having ~26% of HD 80606's surface cov- 
ered with spots. Even if the systematic trends we see in our 
transit light curves are due to spots coming in and out of 
view on the surface of the star, the percent of the stellar 
surface covered by spots is unlikely to be as much as 26%. 
Furthermore, if the star was this spotted, we should also see 
a difference in the flux ratio between the two bluest band- 
passes of over lxlO -4 , yet the difference we observe is less 
than ~6xl0" 5 . 

We conclude that it is possible for spots to account for 
some of the variations we measure, but that HD 80606 is 
very unlikely to be spotted enou gh to cause the mag nitude 
of variations we measure. In fact. I Wright et all |2004h mea- 
sured values of Shk = 0.149 and log R H k — —5.09 for HD 
80606, which indicat e that the star is quite inactive. Also, 
iHebrard et all (|2010T } monitored HD 80606 and determined 
it is not an active star. Specifically, they estimate that the 
star is photometrically stable at the level of a few mmag in 
the optical range on the timescale of several weeks. Despite 
these statements, they still attribute the bump in their light 
curve to a spot on the stellar surface. Considering the pre- 
cision of our observations (much better than 1 mmag), it is 
possible that we observed flux variations that they did not 
have the precision to. 

As a final comment, we note that spots could affect 
the normalization of the overall spectrum, as the spectrum 
could need to be scaled downward (i.e., decrease the flux ra- 
tios or increase the transit depths) to account for the effect 
of spots. However, the shape of the spectrum would remain 
the same, unless over ~26% of the star's surface was covered 
with spots during transit. As noted by the referee, a large, 
long-lived polar spot could exist on HD 80606, which would 
not induce large photometric variations but could still affect 
our photometry. Or, both HD 80606 and HD 80607 or HD 
80607 alone could be spotted and cause the observed varia- 
tions. Due to the possible variable nature of HD 80606 (or 
HD 80607), we encourage future out-of-transit observations 
of HD 80606 and HD 80607 to determine if such variability 
is common. 



5 CONCLUSION 

In summary, our observations do not match existing mod- 
els, due to two basic observations. We find a large change 
in apparent planet radius with wavelength, but do not ob- 
serve a significant difference where the K I line core would 
be expected. Our observations place a strong limit on the 
strength of the line core (unless it has been Doppler shifted 
by > 100 km s _1 ), yet imply large variations in radius over 
wavelengths usually dominated by K I absorption. In the ab- 
sence of other viable absorbers, absorption by K I remains 
the most viable explanation. The atmospheric scale height of 
HD 80606b at transit (~20 km) is significantly smaller than 
that of HD 209458b and HD 189733b, y et the variation i n 
radius is larger than that of HD 209458b (|Sing et alj|2008al ). 
One possible model is absorption by potassium that is part 
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of a high speed wind coming off the exosphere. While high 
speed winds have been observed for other exoplanets, the 
mechanism for powering such winds is unclear. We encour- 
age further theoretical investigations to improve models for 
transmission spectroscopy of exoplanet exospheres in gen- 
eral and the specific challenge of HD 80606b. 

Finally, we have investigated several potential sources 
of systematic effects. There is no simple or obvious source 
causing the systematics in our data. Further, any system- 
atics introduced by the sources we have investigated here 
produce neither the same signature as our observed spec- 
trum, nor the same magnitude of difference as that of our 
measured flux ratios. While we are confident that none of 
these possible sources of systematic effects cause the shape 
of our observed spectrum, we still allow the possibility that 
one or some combination of these systematics may affect 
our measurements and/or the overall normalization of the 
observed spectrum. We also acknowledge that the target was 
observed at a slightly different set of wavelengths in January 
as compared to April. While the difference in wavelengths is 
small (~0.1-nm), there is still the possibility that this could 
result in small differences in either the telluric absorption 
or stellar spectra, which in turn could cause the observed 
spectrum that we have attributed to absorption from the 
atmosphere of HD 80606b. As a final note, we highly en- 
courage follow-up transit observations of HD 80606b to con- 
firm the signal measured here. We note that the next partial 
transit observable from La Palma occurs on 20f2 March 3, 
during which observations pre-transit through the complete 
first half of the transit will be possible. 



6 FUTURE PROSPECTS 

Future transit observations at wavelengths around K I in HD 
80606b are possible, but require considerable patience due 
to the long orbital period (fff days). Observations of the 
transit depth around other absorbing species could test the 
exosphere and wind models. Similar observations of other 
planets would enable a comparison of K I strength in both 
the wing and core as a function of star and planet prop- 
erties. We note that shortly before submission, we became 
aware of i ndependent , but similar, observations of another 
exoplanet (Si ng et al.l 1201 lh . Both these and future obser- 
vations of additional exoplanets will enable comparisons of 
the atmospheric composition and structure, as well as stud- 
ies of potential correlations with other planet or host star 
properties. Such observations would also help improve the 
interpretation of the existing HD 80606b observations. 

Currently, only the OSIRIS red tunable filter (651-934.5 
nm) is available at the GTC. Once the blue tunable filter is 
available for scientific observations, it will be possible to ob- 
serve additional atmospheric features, including the Na I fea- 
ture previously detected for HD 209458b and HD 189733b. 
The large aperture of the GTC makes it practical to perform 
similar observations of several fainter host stars. Thus, we 
look forward to future observations of a large sample of tran- 
siting planets. The striking diversity of exoplanets suggests 
that it will be fruitful to compare Na I and K I observations 
to identify trends with stellar and planet properties. 

Despite the complex interpretation of these observa- 
tions, the very-high-precision obtained with the OSIRIS 



narrow-band tunable filter imager opens up new avenues of 
research for large ground-based observatori es. Indeed, the 
meas ured precision exceeds that of Spitzer (|Hebrard et al.l 
120101 ) and even the Hubble Space Tele scope observations for 
the given bandpass l|Pont et al.ll2008l ). Thus, ground-based 
observations can now characterize the atmospheres of giant 
planets using spectrophotometry. The photometric precision 
is also sufficient to measure emitted and/or scattered light 
during occultation at multiple near-infrared wavelengths 
that could improve constraints on atmosphere models of 
short-period giant planets. By providing high-precision pho- 
tometry at multiple wavelengths during a single transit, the 
technique could also contribute to the confirmation of tran- 
siting _rjlane^_cajididates, such as those identified by Ke- 
pler (|Borucki et al.| [2010). The technique could also improve 
measurements of th e impact parameter and thus orbital 
|Colon fc Fordl l 2009). This would be particul arly valuable 
for sy stems with multiple transiting planets (jSteffen et al.l 
|2010| ) for which the orbital e volution depends on the re la- 
tive inclination of the orbits (|Ragozzine fc Hohn an 2010). 

Since all Neptune and super-Earth-sized planets will 
have relatively low surface gravities, they can make good tar- 
gets for transmission spectroscopy. Despite a smaller transit 
depth than giant planets, the potentially large atmospheric 
scale height can lead to a sub stantial signal in transmis- 
sion l|Charb onncau et all [2009). particularly for Neptune- 
sized planets orbiting sub-solar-mass stars and/or super- 
Earth-sized planets orbiting low-mass stars. Previously, it 
has generally been assumed that the Earth's atmosphere 
will prevent ground-based facilities from achieving the high- 
precisions necessary to measure biomarkers on super-Earth- 
sized planets and that the James Webb Space Telescope will 
provide the first opportunity to cha racterize atmospheres 
of super-Earths (|Deming et al.l 120091 ) . If the challenges of 
Earth's atmosphere could be overcome, then ground-based 
observatories have several advantages (e.g., much larger col- 
lecting area, more modern and sophisticated instrumen- 
tation, ability to adjust and upgrade instruments). These 
observations demonstrate that ground-based narrow-band 
photometry on large telescopes can deliver the precision 
necessary to characterize super-Earth-size planets around 
bright, nearby, small stars. We encourage astronomers to 
consider a future generation of instruments specifically de- 
signed for high-precision transit observations, which may al- 
low the characterization of super-Earth-sized planets in up- 
coming large ground-based observatories [e.g., Giant Magel- 
lan Telescope (GMT), Thirty Meter Telescope (TMT), and 
Extremely Large Telescope (ELT)]. 
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